Search for: All records

Reconstructions of the paleoclimate indicate that ancient climatic fluctuations on Earth are often correlated with variations in its orbital elements. However, the chaos inherent in the solar system’s orbital evolution prevents numerical simulations from confidently predicting Earth’s past orbital evolution beyond 50–100 Myr. Gravitational interactions among the Sun’s planets and asteroids are believed to set this limiting time horizon, but most prior works approximate the solar system as an isolated system and neglect our surrounding Galaxy. Here we present simulations that include the Sun’s nearby stellar population, and we find that close-passing field stars alter our entire planetary system’s orbital evolution via their gravitational perturbations on the giant planets. This shortens the timespan over which Earth’s orbital evolution can be definitively known by a further ∼10%. In particular, in simulations that include an exceptionally close passage of the Sun-like star HD 7977 2.8 Myr ago, new sequences of Earth’s orbital evolution become possible in epochs before ∼50 Myr ago, which includes the Paleocene–Eocene Thermal Maximum. Thus, simulations predicting Earth’s past orbital evolution before ∼50 Myr ago must consider the additional uncertainty from passing stars, which can open new regimes of past orbital evolution not seen in previous modeling efforts.

Given the inexorable increase in the Sun’s luminosity, Earth will exit the habitable zone in ∼1 Gyr. There is a negligible chance that Earth’s orbit will change during that time through internal Solar System dynamics. However, there is a ∼ 1 per cent chance per Gyr that a star will pass within 100 au of the Sun. Here, we use N-body simulations to evaluate the possible evolutionary pathways of the planets under the perturbation from a close stellar passage. We find a ∼ 92 per cent chance that all eight planets will survive on orbits similar to their current ones if a star passes within 100 au of the Sun. Yet a passing star may disrupt the Solar System, by directly perturbing the planets’ orbits or by triggering a dynamical instability. Mercury is the most fragile, with a destruction rate (usually via collision with the Sun) higher than that of the four giant planets combined. The most probable destructive pathways for Earth are to undergo a giant impact (with the Moon or Venus) or to collide with the Sun. Each planet may find itself on a very different orbit than its present-day one, in some cases with high eccentricities or inclinations. There is a small chance that Earth could end up on a more distant (colder) orbit, through re-shuffling of the system’s orbital architecture, ejection into interstellar space (or into the Oort cloud), or capture by the passing star. We quantify plausible outcomes for the post-flyby Solar System.

We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of ∼25 mas (∼8 au). This exquisite resolution enabled us to resolve the compact inner-disk structures surrounding the components of each multiple system and to determine the projected 3D orientation of the disks (position angle and inclination) to high precision. We performed a statistical analysis on the relative alignment of disk pairs to determine whether the disks are preferentially aligned or randomly distributed. We considered three subsamples of the observations selected by the companion separationsa< 100 au,a> 500 au, anda< 10,000 au. We found for the compact (<100 au) subsample, the distribution of orientation angles is best described by an underlying distribution of preferentially aligned sources (within 30°) but does not rule out distributions with 40% misaligned sources. The wide companion (>500 au) subsample appears to be consistent with a distribution of 40%–80% preferentially aligned sources. Similarly, the full sample of systems with companions (a< 10,000 au) is most consistent with a fractional ratio of at most 80% preferentially aligned sources and rules out purely randomly aligned distributions. Thus, our results imply the compact sources (<100 au) and the wide companions (>500 au) are statistically different.

Dynamical instabilities among giant planets are thought to be nearly ubiquitous and culminate in the ejection of one or more planets into interstellar space. Here, we perform N-body simulations of dynamical instabilities while accounting for torques from the galactic tidal field. We find that a fraction of planets that would otherwise have been ejected are instead trapped on very wide orbits analogous to those of Oort cloud comets. The fraction of ejected planets that are trapped ranges from 1 to 10 per cent, depending on the initial planetary mass distribution. The local galactic density has a modest effect on the trapping efficiency and the orbital radii of trapped planets. The majority of Oort cloud planets survive for Gyr time-scales. Taking into account the demographics of exoplanets, we estimate that one in every 200–3000 stars could host an Oort cloud planet. This value is likely an overestimate, as we do not account for instabilities that take place at early enough times to be affected by their host stars’ birth cluster or planet stripping from passing stars. If the Solar system’s dynamical instability happened after birth cluster dissolution, there is a ∼7 per cent chance that an ice giant was captured in the Sun’s Oort cloud.

The majority of binary star systems that host exoplanets will spend the first portion of their lives within a star-forming cluster that may drive dynamical evolution of the binary-planet system. We perform numerical simulations of S-type planets, with masses and orbital architecture analogous to the Solar system’s four gas giants, orbiting within the influence of a $0.5\, \mathrm{M}_{\odot }$ binary companion. The binary-planet system is integrated simultaneously with an embedded stellar cluster environment. ∼10 per cent of our planetary systems are destabilized when perturbations from our cluster environment drive the binary periastron towards the planets. This destabilization occurs despite all of our systems being initialized with binary orbits that would allow stable planets in the absence of the cluster. The planet–planet scattering triggered in our systems typically results in the loss of lower mass planets and the excitement of the eccentricities of surviving higher mass planets. Many of our planetary systems that go unstable also lose their binary companions prior to cluster dispersal and can therefore masquerade as hosts of eccentric exoplanets that have spent their entire histories as isolated stars. The cluster-driven binary orbital evolution in our simulations can also generate planetary systems with misaligned spin–orbit angles. This is typically done as the planetary system precesses as a rigid disc under the influence of an inclined binary, and those systems with the highest spin–orbit angles should often retain their binary companion and possess multiple surviving planets.